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Abstract 

We present a comprehensive study of the massless scalar field perturbation in the Reissner-Nordstrom 
Anti-de Sitter (RNAdS) spacetime and compute its quasinormal modes (QNM). For the lowest lying mode, 
we confirm and extend the dependence of the QNM frequencies on the black hole charge got in previous works. 
In near extreme limit, under the scalar perturbation we find that the imaginary part of the frequency tends 
to zero, which is consistent with the previous conjecture based on electromagnetic and axial gravitational 
perturbations. For the extreme value of the charge, the asymptotic field decay is dominated by a power-law 
tail, which shows that the extreme black hole can still be stable to scalar perturbations. We also study the 
higher overtones for the RNAdS black hole and find large variations of QNM frequencies with the overtone 
number and black hole charge. The nontrivial dependence of frequencies on the angular index £ is also 
discussed. 
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I. INTRODUCTION 



It is well known that under perturbations the surrounding geometry of a black hole will experi- 
ence damped oscillations with frequencies and damping times entirely determined by the black hole 
parameters. These oscillations are called "quasinormal modes" (QNM). It carries unique finger- 
prints of black holes which would lead to the direct identification of the black hole existence. Due 
to the potential astrophysical interest, extensive study of QNM of black holes in asymptotically 
flat spacetimes have been performed for the last few decades (for comprehensive reviews see P, 
and references therein). Considering the case when the black hole is immersed in an expanding 
universe, the QNMs of black holes in de Sitter space have also been investigated 0,0]. 

Motivated by the recent discovery of the Anti-de Sitter/Conformal Field Theory (AdS/CFT) 
correspondence, the investigation of QNMs in Anti-de Sitter (AdS) spacetimes became appealing 
in the past years. It was argued that the QNMs of AdS black holes have direct interpretation in 
terms of the dual conformal field theory. The first study of the QNMs in AdS spaces was performed 
by Chan and Mann Subsequently, Horowitz and Hubeny suggested a numerical method to 
calculate the QN frequencies directly and made a systematic investigation of QNMs for scalar 
perturbation on the background of Schwarzschild AdS (SAdS) black holes Q]. Considering that 
the RNAdS solution provides a better framework than the SAdS geometry and may contribute 
significantly to our understanding of space and time, the Horowitz-Hubeny numerical method 
wa. genera..ed to the study of QNM. of RNAdS black holes . Q and late, crosschecked by 
using the time evolution approach In addition to the scalar perturbation, gravitational and 
electromagnetic perturbations in AdS black holes have also attracted attention |9|,|l0|]. Other works 
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on QNMs in AdS spacetimes can be found in 

Recently in llCl Berti and Kokkotas used the frequency-domain method and restudied the scalar 

PQ 

perturbation in RNAdS black holes. They verified most of our previous numerical results in [3,13]) 
however, disaproved the property of the real part of the QN frequency with the increase of the 
black hole charge. Our first aim in this paper is to check their results by using both improved 
Horowitz-Hubeny method and time evolution approach. 

As was pointed out in 7] and later supported in the Horowitz-Hubeny method breaks 
down for large values of the charge. To study the QNMs in the near extreme and extreme RNAdS 
backgrounds, we need to count on time evolution approach. Employing an improved numerical 
method, we will show that the problem with minor instabilities in the form of "plateaus", which 
were observed in f^, is overcome. We obtain the precise QNMs behavior in the highly charged 
RNAdS black holes. 



In addition to the study of the lowest lying QNMs, it is interesting to study the higher overtone 

n 

QN frequencies for scalar perturbations. The first attempt was carried out in 1^. In the present 
work we would like to extend the study to the RNAdS backgrounds. 

It was argued that the dependence of the QN frequencies on the angular index £ is extremely 
weak .10,]. This was also claimed in [l^ . Using our numerical results we will show that this weak 
dependence on the angular index is not trivial. 

The plan of the paper is as follows. In Sec. II we review the background metric, display the 
scalar wave equations and briefly present our numerical methods. In Sec. Ill we describe in detail 
our numerical results. The conclusions and discussions will be presented in Sec. IV. 



II. EQUATIONS AND NUMERICAL METHODS 
A. Geometry and fields 

The metric describing a charged, asymptotically Anti-de Sitter spherical black hole, written in 
spherical coordinates, is given by 

ds^ = -h{r)dt^ + h{r)-^dr^ + r^{de^ + sin^ Od^'^) , (1) 

where the function h{r) is 



2m ^ Q2 ^^2 



hir) = l- — + ^-—- (2) 



We are assuming a negative cosmological constant, usually written as A = —3/R'^. The integration 
constants m and Q are the black hole mass and electric charge respectively. The extreme value 
of the black hole charge, Qmax, is given by the function of the event horizon radius in the form 
Qlax = rl{l+Srl/R^). 

The spacetime causal structure depends on the zeros of h{r). Changing the parameters m, Q 
and R, the function h{r) may have none, one or two positive zeros. In the Reissner-Nordstrom- 
Anti-de Sitter case, h{r) has two simple real, positive roots (r-|_ and r_) and two complex roots (ri 
and r^, star denoting complex conjugation). In the so-called extreme Reissner-Nordstrom- Anti-de 
Sitter case, h{r) has one double positive zero (also denoted r+) and two complex roots (r^^* and 
(rf*)*). The hor izons r_ and r+ with r_ < r+, are called Cauchy and event horizons respectively. 
In this work we will treat the dynamics of fields in the black hole exterior, the submanifold given 
by the patch r+ = {{t, r, 9, (j)),r > r+}. 
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In the region T-|-, we define a "tortoise coordinate" r*(r) in the usual way, 

dr 



r*{r) 



h{r) 



(3) 



Assuming h{r) has two distinct positive roots (non-extreme case), the tortoise coordinate is given 
by 



r*(r) 



2k, 



In (r — r-i-) In (r — r_) — A ln(r^ + pr + < 



+ 



2{B + Ap) 
yjAq-p"^ 



2k 



arctan 



2r + p \ vr 



y/4:q-p^l 2 



with the constants k+, k_, A, B, p and g given by 

1 R^rl 



2k+ (r+ - r_)(3r^ + + 2r+r_ + R'^) ' 



(4) 



(5) 



2k;_ (r+ - r_)(3r2 + + 2r+r_ + ' 



(6) 



A 



R^{r+ + r_)(r| + ri + 2r+r_ + i?^) 
2{3rl +r'i + 2r+r_ + i?2)(3^2 + ^2 _^ 2r+r_ + ' 



(7) 
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(3ri + ri + 2r+r_ + i?2)(3^2 + ^2 _^ 2r+r_ + R"^) ' 



(8) 



p = r_|_ + r_ , g = i?^ + r^ + ri + r4. . (9) 

The constant of integration in the expression Q was chosen so that hmj._>oo r*(r) = 0. 
Consider now a scalar perturbation field ^ obeying the massless Klein-Gordon equation 

= . (10) 

The usual separation of variables in terms of a radial field and a spherical harmonic Y£ ,„(0, ip), 

'^ = Y,--^it,r)YeUe,cf>) , (11) 



leads to Schrodinger-type equations in the tortoise coordinate for each value of i. Introducing the 
null coordinates u = t — r* and v = t + r*, the field equation is given by 



dudv 



(12) 
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where the effective potential V is 



V{r) = h{r) 



+ 1) 2m 2Q^ 2 



+ ^3 ^4 ' 



(13) 



Wave equation (|12j) is useful to study the time evolution of the scalar perturbation, in the context 
of an initial characteristic value problem, explored in this work. 

In terms of the ingoing Eddington coordinates (v, r) and separating the scalar field in a product 
form as 

^ = Y^^-i;{r)YUO,<P)e-'''\ (14) 

im 

the minimally-coupled scalar wave equation (jlUj) may thereby be reduced to 

/i(r)^^^ + [h'{r) - 2iu] - V{r)ij{r) = , (15) 

where V{r) = V{r)/h{r) = h' {r) /r + l[l + I) /r"^ . 

Introducing x = 1/r, Ea. (|15jl can be re-expressed as Eqs. (15-18) in Q]- These equations are 
appropriate to directly obtain the QN frequencies using the Horowitz-Hubeny method. 



B. Numerical methods 

The quasinormal modes in AdS spacetimes are usually defined as solutions of the relevant wave 
equations characterized by purely ingoing waves at the black hole event horizon and vanishing of 
the perturbation at radial infinity. However, we would mention that in the context of the AdS/CFT 
correspondence, there is no consensus on the "correct"' QNM boundary conditions 
use two different numerical methods to solve the wave equations. 
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. We will 



The first met 



rod is the Horowitz-Hubeny method, which has been used extensively in previous 
papers 0, Q, 0, 3. It was found in 7] that with the increase of black hole charge in RNAdS 
spacetimes, a large number of terms in the partial sum to reduce the relative error in the 
computation of QN frequencies is needed. This takes a lot of computer time and cannot easily be 
performed for a sum of order N > 300. We have to use the trial and error method proposed in [ll| 
to truncate the sum to some large A^. However, with big A^ we adopt, initial tiny error may grow 
through recursion relations. We have to improve the precision of all input data with the help of a 
built-in function of mathematica. In addition high precision of recurrence relation is also required. 
After making these improvements, we can get precise QN frequencies in the lowest lying mode with 
the increase of black hole charge and also high overtones. 
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However, the Horowitz-Hubeny method breaks down for large values of the charge for reasons 
disclosed in j?]. Time evolution approach do not suffer the problem and can be counted on for 
the study of the highly charged RNAdS black holes. The usual implementation of the method for 
two-dimensional d'Alembertians was introduced in and later used in several other contexts 



It was observed in _8] that, for the RNAdS geometries, this usual method develop minor insta- 
bilities with large values of r_|_ in the form of "plateaus". To overcome this problem we have used 
another discretization for equation (|12|) . in addition to the usual one in 



given by 



A2 

1 - —ViS) 
16 ^ ^ 

A2 



lb 



(16) 



The points N, S, W and E are defined as usual: N = {u+A,v + A), W = {u+A,v), E = {u,v+A) 
and S = {u, v). The local truncation error is of the order of O(A^). Although this new discretization 
is more time consuming, we have observed, on a empirical basis, that it is more stable for fast 
decaying fields. 

In presenting the numerical results, we will set the AdS radius R = 1 (always) and the black hole 
radius r_|_ = 100 (for the most part). In the following sections we will first study the properties of 



the lowest QNMs for lowly charged black holes, verifying and extending results found in 
Then we will present our results for lowest lying modes for highly charged holes. We will also 
discuss high overtones of the scalar perturbation for RNAdS black holes and show the nontrivial 
angular index influence. 



III. NUMERICAL RESULTS 

A. Overview of the numerical results 



Employing the Horowitz-Hubeny method to the high precision, we obtained two stable classes 
of solutions for the perturbation frequencies. In the first class, oscillation modes (OM) solutions, 
the frequencies have both the real and imaginary parts: lo{OM) = ujr{OM) + iuji{OM), with 
ujf({OM) 7^ 0. In the second class of solutions, the non-oscillation modes (NOM) solutions, the 
frequencies are purely imaginary: uj{NOM) = iujj{NOM). In both cases, the imaginary part of 
the frequency is always negative, which indicates that the solutions are stable. 

This scenario is consistent with the results obtained with time evolution approach. From the 
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Figure 1: Illustration of the filtering process, (top) Data obtained from the characteristic algorithm, (bot- 
tom) The x^-fitting made in the characteristic data subtracted from the Horowitz-Hubeny NOM solution. 
The results obtained from this fitting are compatible with the Horowitz-Hubeny OM results. In the graphs 
r+ = 100, Q/Qmax = 0.41, R = 1, £ = and n = 0. 



results obtained from the characteristic integration routine, it is possible to estimate with high 
precision the oscillatory and exponential decay parameters using a non-linear fitting based in a 
analysis. Comparing to the results obtained by using the time evolution approach, we found that 
modes with frequencies choosing by min is consistent 

with the results obtained by time evolution approach. We emphasize that the numerical concor- 
dance is excellent. Comparing to the time evolution approach can help us find the criterion to 
determine which class of solutions got by Horowitz-Hubney method dominates the decay of the 
scalar field. We will use this criterion to order the modes. 

On the other hand, no matter the criterion used to order the modes, both classes of solutions 
are physically relevant. The characteristic solutions of the wave equation show the asymptotic 
behavior of the field propagation. We observe two distinct "phases" with the increase of the 
electric charge. For small Q, the scalar field decays exponentially and oscillates. The frequencies 
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Table I: Values for the real and imaginary frequencies for the fundamental mode {n — 0), obtained from 
Horowitz-Hubeny method (OM solutions) and filtered characteristic data, for r+ = 100, R = 1, £ = 0. 



Q / Qmax 


Horowitz- 

LORiOM) 


-Hubeny (OM solutions) 
uJi{OM) 


Frequencies after filtering 
ujR (error) uji (error) 


0.39 


140.74 


-340.47 


138.27 (1.75%) 


-328.98 (3.38%) 


0.40 


144.54 


-344.62 


146.96 (1.67%) 


-355. 71( 3.22%) 


0.41 


147.96 


-347.49 


148.78 ( 0.55%) 


-347.52 ( 0.009%) 


0.42 


150.77 


-349.57 


151.10 ( 0.22%) 


-348.96 ( 0.17%) 


0.43 


152.97 


-351.20 


153.71 ( 0.48%) 


-350.88 ( 0.09%) 


0.44 


154.61 


-352.59 


156.02 ( 0.90%) 


-352.70 ( 0.03%) 



in this case are compatible with the OM Horowitz-Hubeny solutions. For Q greater than some 
critical value Qc, the decay is purely exponential. The exponential coefficients are compatible with 
the NOM-Horowitz-Hubeny solutions. 

But, even for large Q/Qmax, the characteristic approach does not exclude the OM solutions. It 
just says that the NOM solution dominates for large t. This does not invalidate the OM solution, 
because \ujf{NOM)\ < [a;/(OM)[ and the oscillatory solutions could be present in the data, but 
hidden by the non-oscillatory solutions. 

To expose the OM solution in the characteristic data, we will apply a filtering procedure: 
take the wave- function ^(v) obtained from the characteristic integration; subtract from it the 
exponential part {ujj{NOM)) obtained from Horowitz-Hubeny method: 

^{v) — y'^iv)- Cexp [loi{NOM)v] ; (17) 

and finally check if the remaining part is an oscillatory-exponential decaying function, compatible 
with the OM solutions. 

Our results show that both modes are present in the characteristic results, although for large 
values of charge, the NOM solutions dominate the asymptotic limit, as indicated in table Q An 
illustration of the x^-ffitings made is shown in Fig. ^ To simplify notation, we will avoid the 
"(OM)" and "(A^OM)", except where it is not clear from the context. 

As far as we have observed, the qualitative aspects of the field decay are independent of the value 
of the event horizon r^. For completeness, we present in table O some non-oscillatory frequencies 
for other values of r+. 



8 




Table II: Values for the non-oscillatory frequencies u!j{NOM) for several values of r+ and £ = 0, 

Q/Qrnax = 0.35, R=l. 



1"+ 


LUi{n = 0) 


ujiin — 1) 


uji{n — 2) 


uji(n — 3) 


uJi{n = 4) 


UJl{7l ~ 5) 


90 


-229.46 


-663.64 


-1105.3 


-1547.9 


-1990.7 


-2433.5 


100 


-437.11 


-1290.6 


-2155.5 


-3020.5 


-3885.5 


-4750.2 


110 


-805.32 


-2448.5 


-4093.9 


-5738.2 


-7381.9 


-9025.3 


120 


-1486.7 


-4562.5 


-7632.0 


-10698.5 


-13763.5 


-16827.8 



B. Fundamental modes 

1. Black holes with small values of charge 

For the black holes with small values of the charge, it was found that increases with Q, while 
u)R decreases with Q [3, O]. The change of the imaginary part of the frequency was supported, 
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Figure 3: Graph of ujr{OM) with Q/Qmax^ using Horowitz-Hubcny method and characteristic integration 
results. For this graph, r+ — 100, R — 1, I = Q and n = 0. 

while the real part of the frequency was found not with monotonically decreasing behavior [lo| . 
We have two complementary methods and expect to produce a more precise picture of the field's 
evolution. 

For black holes with small values of charge, the scalar field decays exponentially and oscillates. 
For Q greater than a critical value, the decay becomes purely exponential. This point can be 
directly seen in the wave function plotting from Fig|3 We read from Fig|2 that the oscillation 
starts to disappear at Qc = 0.3895Qmax (the first several rings in the plot Qc = 0.3895Qmax are 
due to initial pulse). 

For Q < Qc the scalar perturbation experiences oscillation. The behavior of the real part of 
QN frequencies with the increase of the charge is shown in FigEl where results obtained from 
Horowitz-Hubeny and characteristic integration methods are presented. 

It is apparent that the different numerical methods have very good agreement. We found indeed 
the minimum of w/j^ at a certain value of charge, Qe = 0.366Qmax from the Horowitz-Hubeny 
method and Qe = 0.364:5Qmax from the characteristic method. This result confirms that observed 
in After the minima of ujr, ojr increases with Q a bit, but not as much as described in 
Starting from Qc^ the NOM starts to donimate. The real part of the frequency of the perturbation 
vanishes and will not reach a local maxima at Q = 0A7AQmax as got in |lo| . 

Now let us turn to discuss the imaginary part of the QN frequencies. Our two numerical methods 
give consistent results. From Fig01 we found before Qe where ui r h as minima, monotonically 
increases with charge. This result agrees to that obtained in l7l llCl . However, there is no change 
of sign of the second derivative of u!i{Q/Qmax) as observed in 10]. \iOj\ reaches its maxima at the 
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Figure 4: Graph of uij (OM and NOM) with Q/Qmax, showing that uJi{NOM) tends to zero as Q tends to 
Qmax- In the graph, r+ = 100, R = 1, £ = and n = 0. 



same value of the charge Qe as ujr arrives at its minima. For Q < Qe, our previous statement is 
still valid, "If we perturb a RNAdS black hole with high charge, the surrounding geometry will not 
'ring' as much and long as that of the black hole with small charge." [3, Q]. For Q > Qe, Fig. 4 
tells us that \uji\ decreases with the increase of charge until Q = Qmax- The turning point of luj 
appears at not so big value of Q is different from the rough argument in 



2. Extreme and near extreme limits 



From FigOland Fig|l]we got the picture that for Q > Qc the asymptotic scalar perturbation is 
dominated by the non-oscillatory modes, while keeps decreasing monotonically and smoothly 
up to Q = Qmax- We learnt from Fig|l]that as the black hole becomes extreme, ujj tends to zero, 
which seems to confirm the conjecture raised in section 3.3.2 of [lol |. 

However, from the tail behavior (Fig|^ we observed that in the extreme limit, the decay of 
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Figure 5: Log-log graphs of the scalar field in the event horizon, (left) Extreme case for different values of 
£, showing power-law tails. If £ = 0, ^' cx u^^, but if ^ > 0, ^' oc v^^. (right) Approach to the extreme limit, 
showing transition from exponential decay to power-law decay {£ = 1). For all curves, 7'-|_ = 100 and R — 1. 



scalar perturbation still exists, though changing from the exponential to power-law decay. This 
result strongly suggests that due to scalar perturbations the extreme RNAdS black hole is stable, 
not marginally unstable as worried in by considering electromagnetic and axial gravitational 
perturbations. 

The reason of the smooth change of the perturbation behavior can be understood as follows. 
In the extreme limit, the horizon r+ is a double root of h{r), and therefore the expression for the 
tortoise coordinate is given by 



C 



r — Ta 



Aext ln(r2 + Pextr + Qext) 



+ 



2{Bext + Apext) 
V^Qext - P'ixt 



arctan 



2r + Pext 
4qext - P\ 



,2 

ext , 



(18) 
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with the constants A^xt, B^xt, C, pext and qext defined as 



R^{<or% + brlR^ + R*) ^ ^ 



d2 2 



(22) 



Globally, the structure of the extreme space-time changes drastically from that of the non- 
extreme However, we will argue that, restricted to the exterior of the black hole, the transition 
from near extreme to extreme is, in a sense, smooth. To parametrize the approach to the extreme 
limit, we introduce the dimensionless variable 



r+ — r_ 



(23) 



In terms of 5, the near extreme limit is given by < 5 ^ 1, and the extreme case is 5 = 0. 

Let us analyse the two first terms in the expression for the tortoise function Eq.@. We will 
expand these expressions in b. For this, we note that 

In (r — r+) In (r — r_) = 



2k 



Ik. 



1 i?2 



In 1 + 



(5 + O (<5) 



(24) 



_(5 6r2+i?2 (6r2+i?2)2 
At this point, it is necessary to consider another limit. If r is not too close to the horizon, or 



more precisely, if <C r — then 



1 



2fi;+ 



In (r — r+) 



1 



2k_ 



In (r — r_ 



2R\^ (4r2 + fi2) 
_ In _ 



(6r2 +i?2)^ 



2^2 



6rl + -R2 y r — r+ 



+ 0(5) 



(25) 
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Figure 6: Log- log graphs of the function T^(r*), from the near extreme limit approaching the extreme limit. 
At a certain point, if (5 > 0, there is a transition from power-law decay to exponential decay, indicated for 
each graph. For the curves, r+ = 100, R = 1 and £ = 0, 1 (left and right figures respectively). 



On the other hand, for r sufficiently close to r+, or more specifically if r-\-5 ^ r — ry. 

1 R^r+ 



2k 



In (r — r+) In (r — r_ 



2k. 



S 6rl + 



In (r — r_| 



(26) 



The near extreme approximations for the other terms in the tortoise functions is simpler, one 
can substitute A A^xt, B — B^xt, P Pext and q q^xu with an error of the order of 5. Using 
the previous results, we find that in the near extreme limit, we have 



r*(r) 



rlxtir) + 0{5) 



(27) 



We see that for a certain value of r, the tortoise function becomes dominated by the power-law 
term, characteristic of the extreme function. When very near the horizon, the function r*(r) is 
dominated by a logarithmic term, which diverges when r — r+. 
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For the effective scalar potential, using the previous resuhs, we find that, if <^ r — , 
including the extreme limit where 5 = 0: 



The difference between i = and £ > in expression (|28() comes from the presence of a non-null 
term £{i + l)/r'^ in the potential. This is reflected in the asymptotic behavior of the wave functions, 
as seen in FigEI^left). 

If 6 > 0, the behavior of the effective potential can be different if the observation point is 
sufficiently close to the event horizon. More precisely, if r — r+ <^ (5r-|_, the exponential decay 
dominates, and 



The different potential asymptotic forms, when Q is close to Qmax, depend on whether we are 
initially near to the black hole horizon or not. This difference will disappear in the black hole 
extreme limit {5 0) as illustrated in Fig. 6. As 6 approaches zero, the potential is more and more 
dominated by the power-law phase, and less relevant the exponential tail is. In the extreme limit, 
there is only a power-law tail. This smooth change in the RNAdS potential explains the smooth 
change in the wave functions from the near extreme to extreme limit, illustrated in FiglSf right). 

C. Higher modes 

The first numerical study of the high overtone QN frequencies for scalar perturbation was done 
in SAdS black hole Q]. We have extended their investigation to RNAdS spacetimes and found 
that "switching on the charge" brings interesting physics. For the higher modes, the time evolution 
method is not practical. However the Horowitz-Hub eny method is still efficient, at least for values 
of the charge not too big. 

For the same value of the charge, both real and imaginary part of QN frequencies increases with 
the overtone number n. This result confirms that observed in 15]. 

There are large variations of both uj^ and loj with the charge and the overtone number. Fig|7| 
tells us that with the increase of Q, ujr increases slower while \ujj\ increases faster for bigger 
overtone number n. 

For ujR, combining Fig|31 and Fig|HJ we see that the local minima of the "wiggle" appear at 
smaller values of the charge with the increase of n. The "wiggle" of ujr becomes shallower with 




(28) 



V{r*) oc e' 



(29) 
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Figure 7: Graphs of ijjfj{OM) and uji{OM) versus n, for Q/Qmax = 0.1 (represented by triangles) and 
Q/Qmax = 0.4 (represented by bullets). The continuous lines indicate the linear fits. With the increase of 
Q, lor{OM) increases slower with n, while \bJi{OM)\ increases faster with n. For the curves, r+ = 100, 
R= 1 and ^ = 0,2, 4, 6, 8, 10. 

the increase of n and disappears for n > 5 as exhibited in Fig|Hl Further with the increase of n, we 
observed that u)fi starts to be zero at bigger values of charge. For the imaginary part of the QN 
frequencies, the relation with n and Q are shown in Figs l9ll01 We found that with the increase of 
n, the maxima of appear at bigger values of black hole charge. Moreover we saw from Fie llUI 
that before the maxima, increases quicker with charge for bigger n. These richer physics are 
brought by introducing the additional parameter Q in RNAdS spacetimes. 

In |Q| it was interestingly found that in the large black hole regime the frequencies become 
evenly spaced for high overtone number n. For lowly charged RNAdS black hole, our result 
confirmed their argument for fixed value of the charge. The spacing between frequencies are 



ujn+i - = 129.9 - 225i, 



(n — i- oo when (5 = 0), 



(30) 
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Figure 8: Graphs of ujij{OM) with Q/Qmax, showing the "wiggle behavior". Values of Q/Qmax in graphs 
indicating when non-oscillating solution dominates. For the curves, r+ = 100 and £ = 0. 



which is the same as (12) in (la] by choosing r_|_ = 100 there, and 



Wn+i — = 120.6 — 233i, (n — > oo when Q = O.lbQr 



(31) 



uJn+i — w„, = 111.1 — 256i, (n ^ oo when Q = 0.3Qmax)- (32) 

We found that choosing bigger values of the charge, the real part in the spacing expression becomes 
smaller, while the imaginary part becomes bigger. 

D. Dependence on the angular index £ 



We have so far discussed only the QNMs with 1 = 0. Increasing i, for the lowest 



we obtained the effect of decreasing of toj and increasing of in RNAdS spacetimes 







.jdng mode, 
though 
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Figure 9: Graphs of loi (OM and NOM solutions) with Q/Qmax for several values of n. The continuous lines 
indicate the oscillatory solution, dotted lines indicate non-oscillatory solutions. For the curves, r+ — 100, 
i? = 1 and ^ = 0. 

the dependence is very weak |lo| . In studying the higher modes, it was argued that the asymptotic 
behavior of the spacing of high overtones is independent of the value of I jicj ] . 

Our numerical computations show that although QNM frequencies depend very weakly on I, the 
dependence is not trivial. As was found in the lowest lying mode, lvr increases and \loi\ decreases 
with the increase of L Besides, a quadratic adjustment fits very well the dependence of the modes 
will as seen in Fig JllltT^ This property holds independently of the values of charge and the 
overtone number we calculated. However we observed that lor increases and \oji\ decreases slower 
with I when n increases. 

IV. CONCLUSIONS AND DISCUSSIONS 



We have presented a comprehensive study of scalar perturbations of a Reissner-Nordstrom Anti- 
de Sitter spacetime, and computed its quasinormal modes. We observed the presence of two families 
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Figure 10: Graphs of uji{OM) with Q/Qmax, showing the increase of |lm(w)| with n. For the curves, 
r+ = 100, i? = 1, ^ = and n = 0. 



of solutions, both physically relevant. For Q smaller than a critical value, the asymptotic decay 
is dominated by oscillating modes, while for larger values of Q, the decay is dominated by purely 
exponential modes. Similar situation was observed in the Reissner-Nordstrom de Sitter (RNdS) 
geometry, where A was the critical parameter [4]. 

For the lowest lying QNMs, our results verify (or sometimes disprove) and extend many results 
obtained in P, Q) la 1^- The real part of the QN frequency does have a local minimum at 
Q/Qmax = 0.366, however when Q > 0.3895Qmaa;, the NOM starts to dominate. For this class of 
quasinormal modes, the real part of frequency vanishes. When selecting the mode class (OM or 
NOM) which dominates the decay, we found that the imaginary part of the QN frequency \uji\ has 
maxima at the same values of charge where lor has minima. Before these maxima, increases 
monotonically with charge (OM solution), while after that \ujj\ keeps on decreasing until the black 
hole extreme limit (NOM solution). We observed that under the scalar perturbation the imaginary 
part of the frequency tends to zero, as previously conjectured in the electromagnatic and axial 
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Figure 11: Graphs of (jjfl(OAf) with £, for several values of n. The bullets indicate the numerical results 
and the lines indicate the quadratic fittings. For the curves, r+ = 100, R = 1 and Q/Qmax = 0.1. 

gravitational perturbations In the extreme limit, the asymptotic field decay is dominated by 
a power-law tail, which suggests that this geometry is stable to scalar perturbations. A similar 
field behavior was observed in the RNdS case ^4] in the approach to A = limit. 

We have done an extensive search for higher overtones n of the QNMs of RNAdS black hole 
corresponding to scalar perturbations. We have observed large variations of QN frequencies with 
the overtone number and black hole charge. Both real and imaginary parts of QN frequencies are 
evenly spaced for high overtone number n at fixed Q. The spacings between frequencies varies with 
the black hole charge. The nontrivial quadratic dependence of QN frequencies on angular index i 
has also been addressed. 
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